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ABSTRACT 

Formulations for second and higher harmonic frequency up & down conversions, as well as multi photon processes 
directly assume summability and divisibility of photons. Quantum mechanical (QM) interpretations are completely 
congruent with these assumptions. However, for linear optical phenomena (interference, diffraction, refraction, material 
dispersion, spectral dispersion, etc.), we have a profound dichotomy. Most optical engineers innovate and analyze all 
optical instruments by propagating pure classical electromagnetic (EM) fields using Maxwell’s equations and gives 
only ‘lip-service’ to the concept “indivisible light quanta”. Further, irrespective of linearity or nonlinearity of the 
phenomena, the final results are always registered through some photo-electric or photo-chemical effects. This is 
mathematically well modeled by a quadratic action (energy absorption) relation. Since QM does not preclude 
divisibility or summability of photons in nonlinear & multi-photon effects, it cannot have any foundational reason 
against these same possibilities in linear optical phenomena. It implies that we must carefully revisit the fundamental 
roots behind all light-matter interaction processes and understand the common origin of “graininess” and “discreteness’ 
of light energy. 
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1. INTRODUCTION 

Motivation for this paper derives from the following questions: Does indivisible photon picture give any 
better understanding of the light-matter interaction in nonlinear optical processes? Or is indivisible photon 
based interpretation a physical reality? In this paper, we will investigate photon picture interpretation in the 
context of second order nonlinear optical processes. 

The era of experimental nonlinear optics began with the invention of laser. It set in motion 
experimental nonlinear optics. Experimental studies to investigate the behavior of high intensity light 
beams (optical radiation) in media were carried out with increased peak powers. It enabled the 
advancement of nonlinear optical processes and allowed new observations such as new optical frequency 
changes in nonlinear optical media (second order nonlinear process or Three Wave Mixing) and variation 
of refractive index with optical intensity (self-focusing). It underscored linearity or nonlinearity is a 
property of the medium and not of the optical radiation. Nonlinear behavior is not exhibited when light 
travels in free space, even when the intensity is extremely high. Various applications of light-matter 
interaction or light interacting with light via a medium came to limelight. Nonlinear optical experiments 
have made profound impact on harnessing energy of photon wave packets for various applications 
including optical communications, remote sensing, and spectroscopy. 

1.1. Nonlinear Optics 

Nonlinear optics [1-5] deals with interaction between electric field and atomic field. It deals with the 
effects of non-linear response of a dielectric medium at the atomic level to the electric fields of an intense 
optical radiation. It encompasses the processes of governing energy transfer between optical fields of 


different frequency or direction as intensity dependent scattering and absorption process. It is observable 
when electric field amplitudes of optical radiation become comparable with coulomb fields near atoms 
(~10 6 V/m). Nonlinear contributions become important at high optical intensities. Intense radiation 
produced by a laser propagating through an optical medium produces changes in the spatial and temporal 
distribution of electric charges which result in the displacement of valence electrons from their normal 
orbits. This perturbation creates electric dipole moments whose manifestation is the induced polarization, 
i.e. induced density of electric dipole moment. The electric field of the propagating high intensity optical 
radiation from a laser causes electric polarization in the medium. This polarization propagates together 
with the electromagnetic field in the form of a polarization wave. 

Nonlinear processes generally do not require the assumption of photons. As such, they are modeled by 
treating atoms as nonlinear oscillators without any new quantization process. The nonlinear oscillator 
behavior is described by Lorentz oscillator model for the equation of motion for an atomic electron in an 
electric field E(t). A perturbative solution of the nonlinear oscillator equation is sought. 

In the case of light-matter interaction, light is generated or absorbed via undulation of material dipoles. 
We are exploring whether the indivisible photon picture provide a better understanding of nonlinear optical 
processes. All material media under strong stimulation by EM waves will display both linear and nonlinear 
properties but the strength of the nonlinear property will depend upon the strength of the stimulation. 
When stimulation frequency matches with quantum transition, the linear processes dominate. Response of 
bulk medium without quantized energy levels corresponds to classical nonlinear optics as in the case of 
second-order nonlinear optics or three wave mixing. When the quantum energy levels of the atoms and 
molecules are involved, then the process is quantum mechanical in nature as in Raman scattering and two- 
photon Fluorescence). Our focus in this paper is the classical second order nonlinear optical processes, 
a bulk phenomenon 


2. LIGHT INDUCED POLARIZATION 

The propagating light in a medium induces polarization as described by Maxwell’s equation. The relation 
between P and E is nonlinear. Expanding the P and E relationship using Taylor’s Series, we can write 

P = s 0 - j-E + 2-d-E 2 +4-/ 3) -E 3 +... 0) 

Where e 0 is Permittivity of Free Space, % is the Electric Susceptibility of the medium (=Scalar Constant), 
and E is the Electric Field , [V/m]. The vectors P and E at any position and time are parallel and 
proportional. The first term, e 0 . X-E, in Eq. 1 represents linear dielectric medium. The second and 
subsequent terms collectively represents the nonlinear dielectric medum, the effects of which attain 
significance. The term 2.d.E 2 is the second order nonlinearity where d is the nonlinear optical coefficient 
that is used instead of x (2) for practical reasons. The Polarization density (Dipole moment per unit volume), 
P is given by P = N*e*x, where N = Number of molecules per unit volume, e = electric charge and x = 
electron displacement. 

The second (lowest) order of nonlinear polarization can arise from a ^ (2) nonlinearity which can occur 
only in crystal materials with a non-centrosymmetric crystal structure. The nonlinear polarization then has 
a component which depends quadratically on the electric field of an incident light wave. The nonlinear 
polarization contains frequency components which are not present in the exciting beam(s). Second 
Harmonic Generation (SHG), Optical Parametric Amplification (OPA), Sum Frequency generation (Up- 
Conversion) (SFG), and Optical Parametric Oscillation (OPO). Parametric processes are usually 
polarization-dependent: the nonlinearity itself is polarization-dependent, and at least in cases with a / (2) 
nonlinearity also the phase matching, because such media exhibit birefringence. 

Third-order Nonlinear Polarization allows four-wave mixing processes. The next higher order of 
nonlinear polarization can arise from a ^(3) nonlinearity that give rise to various phenomena: The Ken- 
effect is a nonlinear modification of the refractive index which leads to phenomena such as self-phase 
modulation, cross-phase modulation, Kerr lensing and four-wave mixing. A delayed nonlinear response is 
associated with Raman scattering and Brillouin scattering. 



All linear and nonlinear oscillations takes place under any stimulation by light waves, only the strength 
of the nonlinear oscillations are extremely weak. Notice that there is no threshold limit that bars higher 
order stimulations. The x (n) terms being weak, we need high intensity fields, normally generated by lasers. 

Under any stimulation all linear and nonlinear oscillations takes place. Notice that there is no 
threshold limit that bars higher order stimulations - just a matter of strength of stimulating fields. But the 
X (n) terms being weak, we need high intensity fields normally generated by lasers. 

The electric field of light can be represented by 

E(t) = E 0 • cos(o) • t) (2) 

From Eq. 1, for second order effect, 

Polarization cc E*(t) = Eq • cos 2 (co- 1)= y[Eg + Eq cos 2 (2co- t) (2) 

For three wave mixing, we have two fields, Ei(t) and E 2 (t). Hence, 

E(t) = [Ej • eosfio, • t) + E 2 • cos(co 2 • t)] (4) 

Then, 

E 2 (t) = [Ej • cos(o 1 • t) + E 2 • cos(ca 2 • t)] 2 

= \[E\ + E 2 + E, 2 • cos 2 (2 • <», • t) + E 2 • cos 2 (2 • g> 2 • t) 

+ E,E 2 • eosfio, +co 2 ) - 1 + E t E 2 • cos(co 1 -ro 2 )* t] 


The nonlinear component of the polarization, P NL = 2-d*E 2 contains components 
0, la \ ,la ) 2 , (o) A + co 2 ), (a), -(d 2 )- Alternately, we can write, 

E(t) = E„ • cos(a> t) = Re{S ffl • e‘ tet } = ~(& a - C* + • e la *) 

Writing, 

S ra = S ra (z).e ik » z 


E(t) = Re{S m ■ e‘ i(k “ z ‘ sX) } = ±. {S m . + g; . 


For three wave mixing process, we have two harmonic components namely, 

E = - {& • + s; ■ e ,(k » z crt) } + - {g, . e' l(k " I " 0,) + gl ■ e l(k ”*“ 8 ‘ ) } 

Where, 2 1 1 s 2 x 2 2 
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Writing, 


E 3 = - {S 3 ■ e ‘ ,(k » x ' fflt) + g; • e i(k “ z '® t) } 
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We obtain three coupled equations that can be succinctly written as, 


^ocd-S\e- ,IAk ’'; 
dz i 


i,j,k = 1,2,3 


( 6 ) 

(7) 

( 8 ) 


(9) 


( 10 ) 


( 11 ) 


For second order nonlinear processes, the induced polarization is proportional to 



Z m E 2 (t) = z m ■ E„ • cos 2 (ro • t) 


( 12 ) 


The polarization is material dipole undulation. Expanding the cosine function, 

Z m E\t) = f j (2) • [E \ +E \ • cos(2o) • t)] (13) 


In the above equation, the first term represents the DC polarization and the second term represents the 
sinusoidal oscillation. This sinusoidal undulation, when in phase with the fundamental wave, can convert a 
part of the fundamental energy into the new second harmonic radiation. This is a bulk dipole undulation 
process and the energy conversion efficiency depends upon the physical length over which the phase 
matching condition remains valid. The walk-olf and other associated issues arise when the wave-front 
vector and the Poynting vector cannot remain collinear for a given crystalline medium. Maxwell’s 
equation clearly implies that the light emitting sources have to execute sinusoidal oscillations to be able to 
emit sinusoidal light waves. Materials cannot generate cosine squared radiation. Eq.1-13 has no quantum 
mechanical phenomenon built into them. The DC polarization term is purely a classical bulk phenomenon. 

3. PHASE MATCHING AND WALK-OFF 


Phase matching provides maximum and sustained harmonic conversion over physical length of the 
nonlinear medium. The nonlinear mixing products can be efficiently accumulated over a greater length of 
crystal only if phase matching is achieved. Because of differential variations in the refractive indices, it is 
possible that the phase can be reversed and the non-linear radiation can get converted back into the 
fundamental stimulating radiation. Some nonlinear effects, however, are either automatically phase- 
matched (e.g. self-phase modulation) or do not require phase matching (e.g. Raman scattering). Many 
phase-sensitive nonlinear processes, in particular parametric processes require phase matching to be very 
precise. Essentially, this means ensuring that a proper phase relationship between the interacting waves 
(for optimum nonlinear frequency conversion) is maintained along the propagation direction, so that the 
amplitude contributions from different locations to the product wave are all in phase and hence phase 
mismatch has to be close to zero. 


For three-wave mixing processes, phase matching conditions involve frequency and momentum 
matching relations. The frequency matching condition is given by 


C°Z = (°1 + °^2 


(14) 


The phase or momentum matching condition in terms of wave vectors is given by 

k 3 =k 1 + k 2 (15) 

In the case of optical parametric oscillations, the following cavity resonance conditions in terms of cavity 
length L, speed of light, c, and integers Ni and N 2 have to be satisfied for signal and idler wavelengths. 
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(16) 


Integration of Eq. 1 1 using slowly varying approximation provides field amplitudes of electric field 
amplitudes. The electric field amplitudes are related to phase mismatch, Ak via a Sine function shown 
below. 



where, Ak = k t + k 2 - k 3 ; Phase Mismatch 


(17) 


The term coherence length L c = rc/(Ak), is often used to describe the interaction length when it decreases to 
40% of its peak value. Finally, Manley-Rowe Relation in Eq. 18 for three wavelengths describes energy 
conservation. 



Phase-matching Techniques- A Brief Review 


The common technique for achieving phase matching in nonlinear crystals is birefringent phase matching 
based on crystal birefringence. The birefringent technique to appropriately compensate for phase velocity 
mismatch comes in many variations. In the case of second order nonlinear effect (or three wave mixing), 
three beams are involved namely, pump, signal and idler. Optical parametric oscillation can be represented 
by pump = signal + idler indicating signal (lower wavelength) and idler (higher wavelength) beams are 
generated from the pump beam. In this context, Type I phase matching is given by e = o + o relation 
where pump beam corresponds to extra-ordinary wave (e) and signal and idler beams correspond to 
ordinary wave (o). Similarly, Type II Phase matching is given by e = o + e for the three beams. Besides 
desired signal wavelength and availability of pump beam wavelength, bulk (physical and birefringence) 
properties of the nonlinear medium dictate the type of phase matching to be chosen for efficient phase 
matching process. Critical phase matching is referred to angular alignment of the crystal (or the beam) is 
used to find a phase-matching configuration. Noncritical phase matching refers to all polarization 
directions are along the crystal axes. The angular position is then not a sensitive parameter and the phase 
mismatch is minimized by adjusting the crystal temperature such that the phase velocities of the interacting 
beams are equal. This technique avoids spatial walk-off that will be discussed later. Also, if the wave 
vectors of all involved beams have the same direction it is referred to as collinear phase matching. In the 
case of different directions, it is known as noncollinear phase matching, In any case, the vector sum of 
the generating beams equals the wave vector of the product beam. A special case is achromatic phase 
matching where at least one of the interacting beams is angularly dispersed so that each frequency 
component of the signal is properly phase-matched. In certain cases, phase matching doesn’t occur or even 
if it occurs it may not be efficient in a given crystal direction due to low nonlinear figure of merit. In such 
cases, Quasi-phase matching technique is employed using engineered materials. In this case, the sign (or 
strength) of the nonlinearity varies periodically and obtained by poling techniques. 

In the context of phase matching, group velocity mismatch and bandwidth limitations assumes 
significance. Under phase matching condition the group velocities of the interacting waves are, in general, 
still not matched; there exists a certain group velocity mismatch, which limits the interaction length for 
pulses and (for a given interaction length) the spectral range (called phase-matching bandwidth) in which 
phase matching is achieved. Also, there is only a finite range of beam angles where phase matching works 
- particularly for critical phase matching called as the angular phase-matching bandwidth 

Walk-off 


Walk-off phenomenon assumes significance in second order nonlinear processes limits harmonic 
conversion efficiency [6]. Spatial Walk-off - the phenomenon that the intensity distribution of a beam in an 
anisotropic crystal drifts away from the direction of the wave vector. 

In an isotropic medium, the transverse intensity distribution propagates along the beam axis as defined 
by the medium k vector whereas in anisotropic (and thus birefringent) crystals, the intensity distribution 
drifts away from the k vector direction. The Poynting vector defines the direction of energy transport, 
whereas the k vector is normal to the wavefronts. This phenomenon, called spatial walk-off, ] birefringent 
walk-off or Poynting vector walk-off (not to be confused with temporal walk-ofi), is associated with some 
finite angle p (called walk-off angle) between the Poynting vector and the k vector 


P = ~ 


1 dn e 

n e d6 


(19) 


Phase matching and Walk-off are purely classical phenomenon based on bulk material properties of 
crystals and doesn’t require the use of any quantum mechanical formulations or concepts (use of 
Schrodinger Equation, 2 nd quantization of EM wave fields, energy transitions levels, etc.). 




4. Photon Picture: As described in the standard literature 


In quantum physics, photon is a quantized elementary particle, the quantum of the electromagnetic field 
and the basic "unit" of light (or all forms of electromagnetic radiation). Some times it is described as a 
monochromatic mode of the vacuum with energy hv , where V is the uniquely defined frequency. Photons 
are considered stable Bosons with zero rest mass, zero electric charge and unit spin. Being Bosons, photons 
occupy same space without interfering with each other [7 ]. 

5. Photons and Nonlinear Optics Processes-A discussion 

Often in literature, people invoke quantum mechanics to describe nonlinear optical processes. Figure 1 
illustrates SHG and down conversion cases. However, illustrations in Figure 1 are extended to other three 
wave mixing processes. Nonlinear optical processes do not require the assumption of photons hence are 
often modeled by treating atoms as nonlinear oscillators. 



Figure 1. Quantum mechanical assumptions used in discussing nonlinear optical processes. 


Photon Based Interpretation - Mathematical formalism or physically reality? 


Often people invoke photon concepts to explain nonlinear wave mixing processes. Is the term Three 
Photon Interaction or mixing a physical reality? 

Classical frequency and phase matching equations given in Eq.14 & 15 are often re-written by 
multiplying them through out by (h/27c) as: 





( 20 ) 
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Also, Manley Rowe relation [Eq.18] describing conservation of energy is written in terms of photon 
number conservation based on photon flux densities (photons/s-m 2 ) 


( 22 ) 


d</> \ _ d<f>. 2 _ d<j> 3 

dz dz dz 


The equality of classical relations are maintained while by multiplying appropriate constants both sides and 
hence to create a quantum mechanical interpretation. This is an ad hoc approach to convert purely 
classically derived relations without any physical explanations as to why one needs to make them look like 
quantum mechanical. The experimental validity of Eq.20-22, or any equation, of course, will remain intact 
when they are multiplied throughout by the same constant. But this process does not help us understand the 
nonlinear processes any deeper than what the classical mathematical relations, Eq.14, 15 & 18 have 
provided us. Nonlinear optical processes are classical bulk phenomena. 

Photons being non-interacting Bosons, photons can occupy same space at the same time and photon- 
photon interaction is negligible. How can a stable elementary particle undergo self-interference? Optical 
radiation of terawatt levels can be focused through a pin hole. In nanophotonics, 600 nm wavelength light 
beams have been propagated through glass fiber of ~100 nm diameter which implies photons can be 
“compressed” from a 1mm beam into a lnm beam. All these properties are quantitatively and very 
precisely modeled by classical wave equation and Huygens-FreslnePs diffraction theory (and consequent 
laser beam propagation theory). Introduction of indivisible photon model not only does not bring out any 
better understanding of the physical processes behind nonlinear optics, it creates further confusion in 
distinguishing between real quantum phenomena like Raman Effect and multi-photon atomic processes and 
pure classical bulk crystalline processes. The key point is that without light-matter interaction, and 
mediated by collective dipole undulations of bulk material and stimulated by incident “fundamental 
frequency” of light, we cannot get new optical frequencies. Only material dipoles can create or annihilate 
E.M radiation. Again, photons being Bosons, cannot interact with each other to create new photons. 
Researchers are looking into various predispositions of light-matter interactions to understand true nature of 
photons [8-11]. 


juantum and Bulk Phenomenon 


Quantum theory allows discrete energy level transitions. In traditional nonlinear optical processes 
(frequency up or down conversion using crystals), no quantum mechanically defined energy-band 
transitions or discrete energy level transitions are involved. Continuous energy transfer from primary 
frequency to secondary frequency or pump to signal and idler takes place over long distances in bulk 
medium inside the nonlinear crystal when phase matching conditions are satisfied, x represents bulk 
characteristics. Quasi-phase matching allows energy conversion over increased lengths. Maker fringes 
confirm the phase matching process. As discussed earlier, walk-off is due to dispersion, pump and signal 
or idler energy vector separation. Walk off is a bulk crystal property. During the energy transfer process, 
if the phases are reversed (180° apart), back conversion takes place (signal and idler energies converts back 
into pump) which is easily understood by classical physics. Walk off and Back Conversion are not 
Quantum mechanical processes. Accordingly, frequency summation and momentum vector summation is 
not a justification to invoke bulk nonlinear processes as a quantum mechanical phenomena. These 
equations are representative formalisms with no underlying Quantum Mechanical physics. In contrast, 



Figure 2. Linear and Nonlinear Quantum mechanical processes. Left - linear resonance fluorescence. Right - nonlinear 
(Raman) molecular transitions. 



CONCLUSIONS 


Only light-matter interaction process can generate new frequencies. Bulk crystalline dipole undulation 
via light-matter interaction allows new frequencies to be generated. The processes could be quantum 
mechanical or classical and could be linear or non-linear. We clam that the traditional nonlinear optical 
processes to generate up- or down-converted new optical frequencies are best explained by classical 
formulation that uses bulk material properties. Conversion efficiency depends upon the strength of the 
stimulating field in the presence of effective phase matching or over long coherence lengths. Phase 
matching and walk-off are classical phenomena. Ad hoc multiplications through out of certain classically 
derived equations with suitable constants to make them look like quantum mechanical relations, do not 
bring out any better understanding of the underlying physical light-matter interaction processes. 
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